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Abstract. We prove existence and uniqueness of the global solution to the Cauchy problem 
on a universe fireworks model with finite total mass at the initial state when the ratio of the 
mass surviving the explosion, the probability of the explosion of fragments and the probability 
function of the velocity change of a surviving particle satisfy the corresponding physical con- 
ditions. Although the nonrelativistic Boltzmann-like equation modeling the universe fireworks 
is mathematically easy, this paper leads rather theoretically to an understanding of how to 
construct contractive mappings in a Banach space for the proof of the existence and uniqueness 
by means of methods taken from the famous work by DiPerna & Lions about the Boltzmann 
equation. We also show both the regularity and the time-asymptotic behavior of solution to the 
Cauchy problem. 



1. Introduction 

We are concerned with existence and uniqueness of the global solution to the Cauchy problem 
on a universe fireworks model (hereafter, UFM) without gravity described by the following 
nonrelativistic Boltzmann-like equation [12] 

%+^ = V(t, *, [ r(t, x, • - T(t, x, Of (1.1) 

Ot OX J R 3 

where / = f(t,x,£) is the total mass distribution with respect to time t, space x and velocity £ 
during the explosion of a huge cloud of matter; f\ = f(t,x,£i); 7] = r}(t,x,£) is the mass ratio 
which is a measure of the mass surviving the explosion; T = T(t, x, £) is the explosion probability 
per unit time for a fragment; P = P(— £ • £i) is the probability density so that P(— £ • £i)eZ£i is 
the probability that a surviving particle will change its velocity £ to £i G <i£i during explosion. 
Assume that T, r\ and P satisfy the following physical conditions 

< r) < 1,0 < r < 1,0<P, (1.2) 

I P(-£-&Rl = lforV£GR 3 , (1.3) 

and that the initial total mass distribution /o = fo(x,£) satisfies 

0</o(x,6G^(R. 3 xR 3 ), (1.4) 

i.e., having finite total mass at the initial state. 

In this paper we shall prove existence and uniqueness of the global solution to the Cauchy 
problem on the nonrelativistic Boltzmann-like equation (|1.1|) of the distribution / with finite 
total mass at the initial state when the ratio rj, the probability T and the probability density 
P satisfy the corresponding physical conditions (jl.2p and (jl.3p . We shall also show both the 
regularity and the time-asymptotic behavior of solution to this problem. 
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The nonrelativistic Boltzmann-like equation is mathematically almost trivial since it is a linear 
transport equation with an integrable kernel and bounded multiplicative coefficients. However, 
this paper leads rather theoretically to an understanding of how to construct contractive map- 
pings in a Banach space by means of methods taken from the famous work by DiPerna & Lions 
[3] about the Boltzmann equation. 

The fireworks model is one of the metagalaxy models which are used to give a physical 
explanation of the expansion. Alfven pQ suggested the first fireworks model for the evolution of 
the universe as an alternative cosmology to big bang and described a possible mechanism for this 
model in an original work in 1983, but many details remain unexplained because of the extremely 
complicated physics involved [12]. Some assumptions on the mechanism of the explosion, such as 
that matter and antimatter are mixed, are invalid in a metagalaxy model given by Alfven & Klein 
[2 J for the highest redshift observed. Also, the highest redshift observed cannot be accounted for 
even if it is assumed that there is an upper limit to the velocities that could be reached by one 
explosion of about 0.7c (c is the light velocity) [9]. A special relativistic model was investigated 
by Laurent & Carlqvist [10] but the evolution of the distributions in configuration space was 
not extended to general relativity by use of their approach. Due to these, Widlund [12] gave 
a fireworks model (defined by equation (II. ip ) without prescribing the details of the explosion 
mechanism and developed the model using a fully general relativistic treatment of the fireworks. 
The nonrelativistic fireworks model, which is given by equation (|l.ip . is a special case of the 
nonrelativistic limit to the general relativistic model considered by Widlund [12]. Widlund has 
explained the observed redshifts only from the start of the study of an approximate equation 
of equation (|1.1|) because it is a differential equation and has an exact solution for many weak 
explosions. Our results mentioned above are not only a kind of reasoning to adopt the research 
method of replacing equation (jl.ip by an approximate equation as Widlund has done, but 
additionally make it possible to do some research work about the redshifts directly from the 
start of the solution to equation (jl.ip . 



One of our results mentioned in the previous section can be described as follows. 

Theorem 1. Equation with physical conditions \1.2t) and $1.3\) has a unique nonnegative 

global distributional solution f(t,x,^) through a total mass distribution fo{x,^) satisfying a finite 



total mass condition ji.^l ) at the initial state, and this solution f(t, x, £) belongs to the Banach 



space L^^+oo);^ 1 ^ 3 x R 3 )). 

Theorem [1] shows existence and uniqueness of solution to the Cauchy problem on UFM without 
gravity no matter how long the fireworks era lasts. 

The proof of Theorem [1] is as follows. We first define a mild solution to equation (jl.ip 
as DiPerna & Lions did in [3j. Then, we establish a Banach space and a mapping of this 
space into itself which uniformly decreases distances, and shows that this mapping has a unique 
fixed point which is a unique global mild solution to equation (jl.ip . This is an application of 
the well-known Banach fixed point theorem. Finally, according to the relations between mild 
solution and distributional solution (e.g., [3], [5]), we know that this mild solution is also a global 
distributional solution. 

Definition 1. Let f = f(t, x, £) be a nonnegative function which belongs to Lj oc ((0, T)xR 3 xR 3 ), 
and assume that for almost all (x,£) £ R 3 x R 3 , 



2. Existence and Uniqueness 






and 




(2.2) 
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where hft denotes, for any measurable function h on (0, +00) x R 3 x R 3 , the following restriction 
to characteristics: 

h*(t,x,0 = h(t,x + t£,0- 
if < T < +00, / is called a local-in-time mild solution to the equation 

<-> 

in the time interval [0, T]. IfT = +00, f is called a global mild solution to equation &2.ty) . 

By Definition [H we have the following result: 

Lemma 1. Suppose that f,h<E Lj oc (R x R 3 x R 3 ), and that for almost all x,£ G R 3 , is 
absolutely continuous with respect to t, ti& G L^ oc (R). Then f is a distributional solution to 
equation S2.S\) if and only if f is a mild solution to equation 112.1 



This result is similar to that given by DiPerna & Lions [3j for the Boltzmann equation. To 
show Theorem [TJ we also have to define a mapping as follows. 

Definition 2. A mapping J is defined as follows: for any f G T and t G [0, T], 

J(f)(t, x, = f (x - 14, - I T(a, x + (a- £)/(<7, x + (a - t)£, i)da 

Jo 

+ f f V (a, x + (a- 0^ (a, x + (a - £i)P(-£ ■ &)/(*, x + (a - t)£, ^i)d^da (2.4) 
Jo Jr3 

where T= {/ = f(t,x,£) : f G L°°((0, T); L X (R 3 x R 3 )) for all T G [0,+oo)} . 

Let the Lebesgue measure be used for all the integrals in this paper. Since the Lebesgue 
measure defined in Rx R 3 x R 3 is c-finite, by use of the Fubini-Tonelli Theorem and the physical 
condition {EH}, we know that /, /# / R8 P(-^ ■ [/ R8 P{-£ • G Lj oc (0, +00) 

for almost every (x, ^) G R 3 x R 3 , i.e., the product P(—£,-£,i)f{- • • ) in equation (12.4|) is integrablc 
with respect to ^1. Thus Definition [2] is valid. 

Take To as a definite time the fireworks era lasts. Then, if < To < i, equation (II. lh has 
a local-in-time mild solution in the time interval [0, To], or say, the total mass distribution of 
the universe fireworks is uniquely determined by the initial total mass data during the fireworks 
era. Indeed, by Definition [2] and with the help of the physical conditions (|1.2[) and (|1.3|) . we can 
easily prove that </(/) G T for all / G and that the following inequality holds: 

™ax \\J(f) - J(/i)||z,i(R3 xR 3) <2T max ||/ - /i|| £ i( R 3 xR3) (2.5) 

for all f,h G T . The inequality A2.5H shows that J is a contractive mapping from T into itself 
with the same norm as in C((0,T );L 1 (R 3 x R 3 )). Therefore there exists a unique element 
h% G T such that J (hi) = hi for a.e. (t,x,£) G (0,Tq) x R 3 x R 3 . In order to prove that hi is 
a local- in-time mild solution to equation (II. ip . it is enough to show that hi is nonnegative for 
almost every (t,z,f) G (0,T ) x R 3 x R 3 . 

Let us take another mapping Ji as follows: Ji(f) = max(0, «/(/)) for all / G F + = {/ : 

a.e. 

f G T and / > 0}. Obviously, J- + is a subset of T. Similarly, we can easily show that the 
mapping Ji maps J- + into itself and is uniformly contractive with the same norm as mentioned 
above. Then there exists a unique element fi G J- + such that Ji(/i) = fi for almost every 
(t,x,£) G [0,T ] x R 3 x R 3 . Thus, if fi = J(fi) for almost every (t,x,£) G [0,T ] x R 3 x R 3 , 
/1 is a local-in-time mild solution to equation (jl.ip through /q in the time interval [0, To], 
and fi a = hi. We will below show that J(fi)°=' fi, or equivalently, J(/i)* °= In fact, by 
equation (|2.4|) . we know that J(/i)* is absolutely continuous with respect to i G [0,To] for 
almost every (£,£) G R 3 x R 3 . We may assume without loss of generality that J(fi)#(t,x,£) 
is continuous for all (t, x, £). To prove that J(/i)* = it suffices to prove that J(/i)* > 0. 
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If J(/i)*(ioi ^Oi £o) < for some point (io> £o) G [0, To] x R 3 x R 3 , we can find some value 
h G [0,T ] such that J{fi)*{t, x , £ ) < for all t G [*i,t ], so that f*{t,x ,£o) = for all 
i G [ti,to]- By equation (I2,4p . we can know that 

> J(/i) # (t,x ,£o) > J(fi)*(h,x ,d ) (t >t> h). 

Repeating the above analysis, we can conclude that </(/i)*(0, xo> £o) < 0, i.e., /o(^O)^o) < 0, 
which is a contradiction. Therefore hx = fx = J {fx) > for almost every {t, x,£) G [0,To] x 
R 3 x R 3 . 

How about To > \l It suffices to consider the To = +oo case, i.e., to prove Theorem [TJ In 
order to do this, we have to establish a subset B of such that J is a contractive mapping 
from B into itself with some norm in the time interval [0, +oo). Now, let us take T a as follows: 

F a = \f: \\f\\ a = sup \e- at [ [ \f{t, x, £)\dxd£ ) < A, f € A 

where a > 0, A = f ||/o||li(r3 x r3)- Then {!F a , \\ ■ \\ a ) is a complete Banach space. By Definition [2] 
and with the help of the physical conditions (II. 2p and (11.30 , it is also easy to see that 

IW)||a<4 (2-6) 

\\J{f)-J{h)\\ a <-\\f-h\\ a . (2.7) 
a 

for all f,h G T a . We know from the two inequalities (|2.6p and (|2.7|) that J is a contractive 
mapping of T a into itself with the norm || • [|„ in the time interval [0, +oo) for any given a > 2. 
Thus, choosing a real ao such that oo > 2 and taking B = J- ao , we can have a unique element 
ho G B such that J {ho) = ho in [0, +oo) x R 3 x R 3 . As stated above about the fixed point 
nonnegativity in the < To < I case, we can easily show that ho is nonnegative for almost every 
(t,x,£) G [0, +oo) x R 3 x R 3 . Therefore equation (II. ip with the physical conditions (|1.2j) and 
(ll.3p has a global mild solution through the total mass distribution fo satisfying the finite total 
mass condition (|1.4p at the initial state. Finally, by Lemma [TJ Theorem [1] holds true. 

We can also construct another contractive mapping to prove Theorem [TJ Let us write F&{i) = 
f r*(cr, x, £)da. Then, using the same device as given by DiPerna & Lions [3] for the Boltzmann 
equation, we have 

Lemma 2. / = f{t,x,£) is a mild solution to equation U.l\) if and only if f = f{t,x,£) satisfies 

f*{t,x,0-f*{s,x,Oe-^~ F# ^ 

= f f V #{a,x,Or*{a,x,C 1 )P{-^^)f*{a,x^ 1 )e^ F# ^- F *^d^da 

Js JR3 

for almost all {x, £) G R 3 x R 3 and all < s < t < +oo. 
Definition 3. A mapping J + is defined as follows: V/ G J- + and \/t G [0,T], 

J + (/)(*, x, = e~ J? [ /o(x _ t£ f £) 

+ f Q + {a,x + {a-t)Z,Zx)e- f ° r{T ' x+{T - m ' 0dT do-}, (2.8) 
Jo 

where 

T+ = {f = f{t,x,0 : / > and f G L°°((0, T); T 1 (R 3 x R 3 )) /or all T G [0,+oo)} , 

q+c*, x, o = / x, e)r(i, x, eopR • x, eo^i- 

7R3 
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Using Lemma [2] and Definition [3l we can give a different and brief proof of Theorem [TJ Let 
us take as follows: 



K = \f-\\f\\a= ™p \e~ at \f(t,x,0\dxdt\<A,fe^ 

where a > 0, A = all/ollL^RSxR 3 )- Then , \\ ■ || a ) is a complete metric space with the distance 
dist(-,-) = || • — • || a . Under the physical conditions (|l-2j) and (II. 3p . it can be easily seen from 
Definition [3] that 

\\J + (f)\\a < A, (2.9) 

ll^-C/) - J^WIU < -11/ - ^IU- (2-io) 

a 

for all f,h£ . By (I2.9P and (|2.10p , it follows that J + is a contractive mapping from into 
itself with the distance || • — • || a in the time interval [0,+oo) for any given a > 1. This implies 
that if we put a > 1, then there exists a unique element go G such that J + (go) = go- Hence, 
by Lemma El Theorem [T] holds true. 

Similarly, we can further show the regularity of solution to the Cauchy problem, that is, 

Theorem 2. Let k be any fixed natural and f a unique nonnegative global distributional solution 
under assumptions in Theorem{J\ Support that r\ and T are in C fc ([0, +oo) x R 3 x R 3 ) and that 
P € C k (R). If f e C k (R 3 x R 3 ), then f G C k ([0,+oo) x R 3 x R 3 ). 

In fact, under all the assumptions in Theorem EJ we can find that J + defined by (12. 8D is still 
a contractive mapping from Pi C k ([0, +oo) x R 3 x R 3 ) into itself for any fixed a > 1, and 
then easily know that Theorem [2] holds. 

Remark 1. In the relativistic case, the total rest mass distribution of UFM without gravity is 
uniquely determined by the initial total rest mass data as in the nonrelativistic case. Since the 
velocity £ in the relativistic case is not greater than the light velocity c (for convenience, write 
c = 1), we have to choose a transformation: £ = = where p is a variable of momentum for 

the relativistic UFM, and regard the total rest mass distribution as that of the three variables 
of time t, space x and momentum p. Thus we can easily show the results mentioned above as 
done in the nonrelativistic case, although the form of the transport operator ^ + —j==J^ of the 

relativistic Boltzmann-like equation [T^] is different from that of the nonrelativistic Boltzmann- 
like equation It is here worthwhile to mention that many results about the Cauchy problems 

relative to the relativistic Boltzmann or Enskog equation have been given in a number of papers 

(e.g., a, m, m,w- 

3. Asymptotic Behavior 

In this section we shall study the time-asymptotic behavior of the solution to the Cauchy 
problem on UFM. We can find that this solution is time-asymptotically convergent to the free 
motion in Z^-norm under some suitable assumptions. 

To do so, we first have to show the following lemma: 

Lemma 3. Let 5 be a positive constant less than 1. In addition to all the assumptions in 
Theorem^, if n and P satisfy the following inequality: 

[ v (t,x,0P(-^i)d^<S (3.1) 

for all the three variables (t, x, then the nonnegative solution f to the Boltzmann-like equation 
i fZHP satisfies 

iff f(t,x t £)dxdZ + (l-8) ff T(t,x,Of(t,x,t)dxdt<0. (3.2) 
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Furthermore, we have 

// f(s,x,£)dxd£ 

J JR 3 xR 3 

<(6-l) f ff T(t, x, Of(t, x, Odtdxd* +11 f (x, Odxdt (3.3) 
for any s 6 (0, +oo). 

Proof. Note that /#(£,£,£) is absolutely continuous with respect to the time t and that the 
solution / satisfies 

§- t f*(t,x,0 = -r#(t,x,0f # (t,x,0 

+v*(t,x,z) I r^x + ^djFK-ei/l^ + tC^i. (3.4) 

J-R? 

Integrating equation (|3.4p over the two variables (x,£) and using assumptions (jl.2p and (|3.ip . 
we can know that (|3.2p holds true. Furthermore, with the help of (jl.4p . integrating (|3.2p over 
the time t from zero to s gives (|3.3p . □ 



Remark 2. T/ie inequality \S.13\i implies that the total mass in UFM does not increase as time 
increases. 

Then, using Lemma EJ we can show that 

Theorem 3. Under all the assumptions in Lemma\^ the solution f(t, x, £) to the Boltzmann-like 
equation hl.l\) converges in L X (R 3 x R 3 ) to foo{t,x,£) as t — > +oo, where 

r+oo 

/oo(t, x, = h(x - - / I>, x + (a - t)£, £)f (a, x + (a- t)£, £)da 

Jo 

r+oo r 

+ / V (a,x + (a-t)^OP(-^i)r(a,x + (a-t)^^)f(a,x + (a-t)0^)dad^ (3.5) 

Jo JR 3 

for any point (t,x,£). 

Proof. Note that / is a nonnegative solution to the Boltzmann-like equation (jl.ip . Then, by 
(|3,5p . we can easily know that 

r r r+oo /• r 

// \f(t,x,0-foo{t,x,0\dxd?< / // r(*,x,Of(<T,x,£)dxd£da 

JjR 3 xR :i Jt iiR 3 xR :i 



+ 



+oo 



T](o-,X,OP(-ttlWl 

R 3 



r(o-,x,€i)f(a,x,£i)dxd£da. (3.6) 



/R 3 xR 3 

By (|3.1|) . it then follows that 

r r r+oo r r 

// \f(t,x,0-fo O (t,x,0\dxdt<(l + 5) / // r(a,x,0/(ff,s,0<ted^- (3.7) 

J^R 3 xR 3 ViR 3 xR 3 

Since the inequality (13. 3p in Lemma [3] implies that the integral on the right side of (13. 7p is 
convergent to zero as t —* +oo, the left side of (|3.7p converges to zero as time goes to infinity. 
This completes our proof of this theorem. □ 

Remark 3. Since foo(t,x,£) defined by \3.5\) describes a free motion in UFM, Theorem^ shows 
that the solution to the Cauchy problem on UFM is time- asymptotically convergent to the free 
motion in L l -norm under some assumptions. It is worth mentioning that this time- asymptotic 
convergence of solution still holds true for the Boltzmann or Enskog equation under some suitable 
assumptions (e.g., [TT] ). 
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